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I. INTRODUCTION 



Understanding equilibrium and noncquilibrium phenomena associated with phase transitions has become 
more and more important in various fields in physics, such as early-time universe, ultra-relativistic heavy-ion 
collisions, ultra-cold atoms, and so on [1, 2]. In a second-order phase transition, characteristic long-range 
fluctuations appear in the order parameter field, and for a quantitative study of static and dynamic critical 
phenomena [3, 4], one needs a field theoretical method which can describe both static and dynamical processes 
involving strong fluctuations. Emergence of long-range fluctuations makes naive perturbation theory break 
down and requires some sort of resummation, such as the method of the renormalization group or the two- 
particle-irreducible (2PI) effective action [5, 6]. 

Recently, the method of the 2PI effective action has received much attention since it can be applied to 
the phenomena in and out of equilibrium on an equal footing [7-9]. In this method, all the self-energy 
contributions for the two-point correlation function are first summed up and then the perturbative expansion 
is carried out in terms of the full two-point correlation function. This is in contrast to the standard method 
of the one-particle-irreduciblc (1PI) effective action in which the perturbative expansion of the diagrams is 
done in terms of the free two-point correlation function. The method of the 2PI effective action systematically 
resums higher order terms in powers of coupling constants, so that it is expected to take into account efficiently 
the large fluctuations near the critical point. 

In the present paper we restrict ourselves to static critical phenomena and leave dynamic critical phenomena 
for future study. As is well known, the most prominent feature of static critical phenomena is universality. 
In other words, several critical exponents which characterize the singularities in the vicinity of the critical 
point are solely determined by symmetry of the system, irrespective of microscopic details. In fact, only two 
of them are independent, and we take n and v to be studied in this paper. They can be read off from the 
two-point correlation function G(x, 0) = (<p(x)<p(0)) of the order parameter field as 

G(x,0)~e-Wt, Z~\T-T C \-" (T>T C ), (2) 

where d is the number of space dimensions (now d = 3), £ is the correlation length, and T c is the critical 
temperature. Namely, rj and v govern the behavior of two-point functions on and off the critical point, 
respectively. In the momentum space, the Fourier transform of Eq. (1) gives the scaling form G(k) ~ |fc| I ' _2 . 

Recently, Afford, Berges, and Cheyne employed the 1/N expansion of the 2PI effective action to compute 
the exponent n of an 0(iV)-symmetric ip 4 theory in three dimensions [7]. They solved the 2PI Schwinger-Dyson 
equation (Kadanoff-Baym equation) [10-12] at the critical point, substituting the scaling form to G(k). It 
was shown that at the next-to-leading order (NLO) in the 1/N expansion this approach remedies the spurious 
divergence of 77 at small N, which is seen in the 1PI 1/N expansion, and leads to an improved estimate already 
for moderate values of N. This success strongly motivated us to compute another exponent v within the 2PI 
effective action. The exponent v is also associated with the critical behavior of the two-point functions. 

It is not straightforward, however, to apply this method to the evaluation of v. Given two nonzero param- 
eters, p and T — T c , one needs to fix the form of G(p; T) in the scaling region, which introduces a technical 
complication to the problem. Fortunately, we notice that the exponent v can be determined from the three- 
point vertex function r^- 2,1 '(x,y;z) ~ (ip(x)(p(y)ip 2 (z)) with two elementary fields, ip, and one composite 
operator, ip 2 evaluated at T = T c [13]. In fact, its Fourier transform (see Eq. (6) for definition) at the critical 
point behaves as 

r^f^^feVlfel 2 -^ (t = t c ). (3) 



v 2 T 

Therefore, if one finds the equation for r^ 2,1 ' in the 2PI formalism, one should be able to compute v at the 
critical point, similarly to the case of rj. 

In this paper we develop the 2PI formalism for the three-point vertex function r^ 2,1 ', and apply the 1/N 
expansion to compute the exponent v. We calculate r' 2 ' 1 ) at the next-to-leading order in the 2PI 1/N 
expansion assuming the scaling form of the correlation function at the critical point. We then extract the 
exponent v according to Eq. (3) and examine whether an improvement similar to the calculation of i] is 
achieved. 
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Computation of the critical exponents has been challenged since 70's, in the e-expansion [14] and 1/N- 
cxpansion [15-18] appoarches in the 1PI effective action formalism. Furthermore, the four-particle-irreducible 
(4PI) effective action has also been applied in [19] to get the higher order terms in the 1/N expansion. These 
methods are utilized to obtain a strict 1/N expansion series eventually. In contrast, our motivation here is to 
examine a possible improvement due to the self-consistent approximation provided in the 2PI formalism. 

This paper is organized as follows. In the next section, we first define our model and then explain how r' 2,1 * 1 
is related to the critical exponent v. The formalism with the 2PI effective action is introduced in the third 
section, where we also derive a self-consistent equation for T 1 - 2 - 1 ' . Then, in the fourth section, we calculate v 
in the 2PI effective action and compare it with the 1PI result, where some complications in the calculations 
are deferred to Appendix. The final section is devoted to a summary of our results and discussions. 



II. THREE-POINT VERTEX FUNCTION r (2,1) AND CRITICAL EXPONENT v 

We consider an O(N) symmetric ip 4 model (ip a = tpi, ■ ■ ■ , cpx) in the three-dimensional Euclidean space. 
The action is given by 



S[<p] = / d 



f->> 



^diip a (x)d t ip a (x) + -^(<Pa(x)ip a {x)) 2 + ^ t ip a (x) if a (x) 



(4) 



where t can be identified as either the mass squared or T — T c , up to renormalization. Roughly speaking, 
the ground state is in a symmetric (broken) phase when t > (< 0), and the transition at t = is of the 
second order. Although the exponents are symmetrical about t = 0, we compute the critical exponents by 
approaching the critical point from the symmetric phase (t > 0) because the vanishing expectation value 
<fi = {(p) = makes the computation technically less involved. 

The two-point correlation function, G a b, and its Fourier transform, G a b, are defined as 

Gab{x,y) = (ip a {x)ip b (y)) 

= Jwf e ^ v)dab{kh (5) 

where translational invariance in the equilibrium state is assumed. Being in the symmetric phase, we treat 
G a b as diagonal and we write G a b = GS a b unless otherwise stated. Similarly the three-point vertex function 
with two elementary fields and one composite field, T^ 2 ' 1 ', and its Fourier transform are defined as 

^ I (0(^ cifc( ^ )c " ip(z " y)F - 1) ( fc '^ fc ^) ■ ( 6 ) 

where summation over the indices a', b' and c is implied and we have assumed translational invariance of an 
equilibrium state. 

Notice that there is a relationship between the two-point function G and the three-point vertex function 
p( 2 ' 1 ). If one regards t as the external field coupled to ^ip 2 , then the differentiation of G with respect to t(z) 
gives [13, 20] 

5G ab (x,y) t _ ,_ v _ ! 2 



St(z) 



= -(M*)<Pb(v)h<f c (*)} • (7) 



Using SG 1 / St = — G 1 (8G/8t)G 1 which follows from the identity G" 1 G = 1, one finds that the definition 
of r' 2 ' 1 ) yields 
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The corresponding equation holds in the momentum space. In particular, in the zero momentum limit, one 
has 

f(V) (0 ,0;0) = ^^). (9) 

At the critical point, t = 0, the exponent i] is determined from the low-momentum behavior of the two-point 
correlation function G(k) ~ |fc| _2+ '', while the exponent v can be obtained from r' 2,1 '(fc,p — k;p) as shown 
in Eq. (3) [15]. This can be explained with the help of the scaling hypothesis applied to T^ 2 - 1 \k^p — k;p). 
Near the critical point, the susceptibility, Xi behaves as G(0) — \ ~ t~~< with the critical exponent 7, which 
immediately implies that 

r< 2 ' 1 >(0,0;0)~t' r - 1 . (10) 
In the scaling hypothesis we assume the existence of a function / and a constant y such that (k = \k\) 



2' 2 

When t 7^ 0, the limit k — > is regular and so is /(0), which yields 

f^' 1 \0,0;0)r^f(Q)C u ^t- vv , (12) 
where the use has been made of £ = t~ v (see Eq. (2)). Comparing Eqs. (10) and (12), one finds 

y=^. (13) 

As we approach the critical point t — > 0, the correlation length ^ diverges, while T^ 2 ' 1 \k/2, k/2; k) is still 
finite as long as k is kept nonzero. Therefore, we must have f(k£) ~ (k£)~ v to find the scaling form at the 
critical point: 

r^Q.^fe)-^-^. (u) 

This is equivalent to Eq. (3) with the aid of the scaling law 7 = v{2 — rj). Diagramatically, it is shown as 




(15) 

where a blob, a wiggly line and a simple line represent r^ 2,1 \ tp 2 and G, respectively. A slash on a simple line 
indicates the amputation. 

III. 2PI EFFECTIVE ACTION 

We give here a minimal review on the 2PI effective action, together with the 1PI effective action for 
comparison. The generating functional or the partition function Z[J] with an external field J is 



Z[J} = J Vipcxp -S[ip}+ J d 3 xJ a (x)<p a (x) 



e" ff M , (16) 



where W[J] is the generating functional for the connected Green's functions. The averaged field is given by 

<f, a ( x ) = { Va ( x )) = ^W. (17) 
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The 1PI effective action Tipi as a function of <f> is obtained by the Legendre transformation of the generating 
functional W[J], 

riPiM = w[j] - f d 3 xj a (x)^\ . (is) 

J oJ a [x) 

Diagramatically, Tipi consists of the vacuum diagrams written in terms of the lines representing the free 
two-point function Go(4>) in the presence of the classical field, 4>. Each of 1PI diagrams does not split into 
two by cutting only one line. 

The ground state is determined by the condition: 6Txpi[</>]/6(j) a (x) = J a (x) = 0, which has a useful form 
for variational analysis. 

In order to obtain the 2PI effective action, we introduce two external fields J and K and define the generating 
functional Z[J,K] as 



Z{J,K} = J Vipcxp -%>] + J d 3 x,J a (x)tp a {x) + J d 3 xd 3 ytp a {x)K ab (x,y)tp b (y) 



e -W[J,K]_ (19) 



Here the generating functional W[J, K] is defined by the last equality. The averaged field (f> a (x) = ((p a (x)) 
and the full propagator (or the correlation function) G ab (x,y) — (<p a {x)ip b (y)) connected are respectively given 

by 

oJa(x) dRab(x,y) 2 

Performing the Legendre transformation of W[J, K] with respect to J and K, we obtain the 2PI effective 
action T2Pi[(j>, G] as a function of <f> and G, 

r 2PI [0,G] = W[J,K] - / d 3 xJ a (x) 6 ^ [J ;* ] [ d 3 xd 3 yK ba {y,x) ^ [ f K] (21) 
J SJ a (x) J SK ab (x,y) 

One can explicitly extract the one-loop contributions from the 2PI effective action in the same manner as 
in the 1PI effective action (but now using the full propagator), yielding the most general and useful form (for 
derivation, see Ref. [9]): 

r 2PI [0,G] = S[<l>] + ^TrlnG- 1 + ^TrG^G + T 2 [0, G] , (22) 

where Tr should be understood as integration over the space coordinates and summation over the field com- 
ponents. The last term T 2 represents contributions from 2PI vacuum diagrams in terms of the full propagator 
G, not of the free propagator Go- 

The ground state is determined by the stationary conditions with respect to <p and G at vanishing external 
fields J = K = 0, and turns out to be the same as in the 1PI effective action, as it should. 

A. Self-consistent equation for G: Kadanoff-Baym equation 

Performing the functional derivative of Eq. (22) with respect to G and setting K = 0, we obtain 



= -lGat(x,v) + ^G-Ux,y 



sr 2 [<j>,G] 



2 ab v »»/ 2 o,a 6 v SG ba (y,x) 

Comparing this with the Schwinger- Dyson equation, G _1 = Gq 1 — S with the proper self-energy E, we find 
that 

X ab { ( f>,G(x,y)} = -2 ^ 2[ ff} (23) 
dG ba (y,x) 



G 



Namely, the functional derivative of T2 is identified with the proper self-energy, which must be 1PI, and 
therefore T2 is 2PI in terms of the full propagator G, as we mentioned above. Thus, we arrive at a self- 
consistent equation for G, the Kadanoff-Baym (KB) equation [10-12]: 

G- b \x, y) = G^ ab (x, y) - £ ab [0, G(x, y)] . (24) 

We remark here the following: if one eliminates G in favor of (f> from T2Pi[4>, G] using Eq. (24) to obtain 
r2Pi[</>, G(4>)] as the functional of </>, one should formally recover the 1PI effective action Tipi^], and therefore 
the ground states in both approaches must be the same. In practice, however, these effective actions are 
different in approximation level because resummation has been done in T2Pi[4>, G(0)] [21]. Introduction of the 
full propagator G satisfying the self-consistent equation (24) provides us of a way to reorganize the expansion 
series in a perturbation theory. 



B. Self-consistent equation for r' 2 ' 1 ' 

Now one can derive the self-consistent equation for I^ 2,1 ) from the KB equation (24) for G. By differentiating 
Eq. (24) with respect to t(x) and using the relation (8), one obtains 

t^(2,i)/ \ r t \ 5^ab[G{x,y)] 

Kb '(x,y;z) = r ,ab{x,y,z) ^— , (25) 

where Tq^^x, y; z) = S(x — z)S(y — z)5 a b- Because we are in the symmetric phase <j> = (cp) = 0, we have 

Gab = G(5 a t,, S a t, = £<5 a &, and rjj? = r' 2 ' 1 ^;,, and thus we deal with the scalar functions without indices, 
hereafter. Applying the chain rule for S = E[0 = 0, G], we can rewrite the second term as 

SE(x,y) f 3 3 SGix^yJ 8T,(x,y) 



d x\d y\- 



6t(z) J ai St(z) 5G(xi, Vl ) 

= - / d 3 x 1 d 3 y 1 d 3 x'd 3 y'G{x u x') ) G(y , Vi) xru \ 

J ot(z) dG{xi,y 1 ) 

= - J d 3 x'd 3 y'D(x, y; x\ y')T^ (x', y';z), (26) 

where we have used -^G = — G (^G _1 ) G, similar to the one used for Eq. (8), and defined the kernel D as 

D(x,y;x',y') = f d 3 x l d 3 y 1 G{x 1 ,x l ) ^ v) G(y', Vl ) . (27) 
J 5G{x 1 ,y 1 ) 

Thus, we write the self-consistent equation for I^ 2 ' 1 ) as 

(x, y; z) = T (x, y;z) + J d 3 x'd 3 y'D(x, y; x', y')^ 2 ' 1 V, 2/'; *) • (28) 
In the momentum space, we have 

r (2ll) (p,g;p + g) = i + J ^ ^ 3 D(p 7q -,p', q ')r^\p', q >; P > + q ') , (29) 

where the Fourier transform of the kernel is defined by 

D(p, <?; p\ q ') = j d 3 xd 3 yd 3 x'd 3 y'c l (- x p - y q+x ' p ' +y ' q ' ) D{x 7 y; x' , y') . (30) 

The kernel D(p, <j; p', q') contains the momentum conservation condition (2ir) 3 5(p + q — p' — </'). We evaluate 
the equation (29) to calculate the vertex function I^ 2 - 1 ) with a given kernel D at the critical point t = 0, and 
determine the exponent v. 
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IV. CRITICAL EXPONENTS FROM 2PI EFFECTIVE ACTION 



As we explained in the previous sections, the exponents rj and v are respectively associated with the two 
point function G and the three-point vertex function I^ 2 ' 1 ) at the critical point. To determine the long distance 
behavior of these functions, one needs to solve the self-consistent equations (i.e., Eqs. (24) and (29)), both of 
which are derived from the 2PI effective action T^. 

Let us first evaluate ^[G] (22) up to the NLO accuracy in the 1/N expansion. The LO and NLO contri- 
butions are respectively 



r LO 

1 2 



7 NLO 



[G] - OO = J d 3 *G aa ( X , X 



[G] 




— Tr In 

2 



)G bb (x,x) , 



S(x - y) + —G ab (x, y)G ab {x, y) 



(31) 



(32) 



where a gray blob indicates a vertex X/N and a line corresponds to G a b = G5 a b- Summation over the repeated 
indices a, b, etc. should be understood. Then a closed loop gives the number of the field components N, and 
thus the first diagram amounts to 0(N 2 /N) = O(N), while the second 0(N°). One obtains the self-energy 
at the NLO by cutting one propagator in these diagrams, which yields in the momentum space 



Sab(p) 



o 



d 3 q 
W) 



3 G{q)5 ab 



d 3 q 



(27T) 



3 Hq)G{p - q)Sab , 



(33) 



where the dashed line corresponds to the sum of bubble chain diagrams, which is denoted by /: 



X + + )o< >oC 



A 



3^V i + | n(p) 
with the one-loop polarization function 

n( P ) = 



d 3 q 



G(p-q)G{q). 



(34) 



(35) 



One should keep in mind that the vertex / is an 0(1/N) quantity. With the NLO self-energy (33), the KB 
equation for G^^p) = G~ 1 {p)8 a b reads 



G-\p) = G^(p) + - 



d : '' ( 



6 J (2tt) 3 



G(q) 



d 3 q 

(27T) 2 



I(q)G(p-q) 



(36) 



A. Calculation of rj 



In Rcf. [7] the critical exponent rj was calculated from the KB equation (36) at the critical point. Let us 
briefly review here how to obtain rj, which is also necessary for the calculation of v. 

Recall that G(p) should become massless at the critical point: G _1 (0) = 0. We can make this condition 
explicit for the KB equation (36) by subtracting the corresponding equation evaluated at p = 0. Then, the 
KB equation at the critical point reads: 

d 3 q 



G-\p) = p 2 



(2tt)s 



HQ) G[p-q)-G{q) 



(37) 
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As one approaches the critical point, the polarization II dominates in the denominator in the scaling region, 
and therefore we can ignore "1" in the denominator of I(p) on the right-hand side of Eq. (34): 

j(p) ~ fn-^p). (38) 

Indeed, in order to investigate the asymptotic behavior of Eq. (37) in the small momentum region, we can 
use the scaling form for G(p) 



with p = \p\ and a cutoff scale A, and find that H(p) is infra-red singular as long as r\ < 1/2: 

U(p)=A( V )^^) 2 \ (40) 

where 

Mr,) - 1 r(j-„)[rm] 2 

After performing the remaining integral with the use of the scaling form (39) and rescaling with A of the 
scaling region, the KB equation (37) reduces to 

p 2 -" = p 2 + B( V )p 2 -" + F> 2 ) , (42) 

where 

= 4/7(1 - 2r]) cos(r/7r) 
U> (3-77)(2-? ? )sin 2 (7?7r/2)A 

and 



" ~ 6^(7,) A P +U[P 



In Eq. (42), terms with p 2 v are dominant at small momentum, p ~ 0, and determine the long distance 
behavior. Equating the coefficients of p 2_?? , one observes that 77 has to satisfy 

1 = B{rf) . (43) 

This gives the NLO result of 77 in the 2PI 1 /N expansion. In Fig. 1 we show the exponent 77 fixed by Eq. (43) 
as a function of N in a solid line, and compare it with the 1PI result written in a dashed line. It is evident 
that the divergence of 77 at N = 0, which is seen in the 1PI result, is now resolved in the 2PI result, and that 
the 2PI result is closer to the experimental values [22]. 

B. Calculation of v 

1. Self- consistent equation for r' 2,1 ^ and iterative solution 

Now, we proceed to the calculation of v in the 2PI formalism up to the NLO in the 1/N expansion. As we 
explained in Sec. II, we will compute v from the three-point vertex function r' 2 ' 1 ^(fc/2, fc/2; k), which satisfies 
the self-consistent equation (29). Thus, the first thing to do is to determine the NLO form of the kernel D in 
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FIG. 1: The critical exponent rj as a function of TV. The 2PI (1PI) result is shown in a solid (dashed) curve. 



Eq. (29). According to the definition of D in Eq. (27), it is obtained by differentiation of the self-energy £ 
with respect to G in the NLO approximation and is written explicitly as 



D(p,q;p\q') = ,<^V + pi 



'•I li* 



-- 6 - Y (P'-P^"J W y 
x G(p')G(q')(27r) 3 S(p + q-p' - q') 

=d + ± Dl . 



I S^G(k - p)G(k - p')I(k)I(k -p-q) 



(44) 



The first term is of 0(1), which is obtained from the first diagram in Eq. (33) by cutting the loop. The second 
and third diagrams are of 0(1/N) (recall that / is of 0(1/N)), which are obtained from the second diagram 
in Eq. (33) by cutting the solid and dashed lines in the loop, respectively. We introduced by the last equality 
the shorthand notations Dq and jjDi , respectively, for 0(1) and 0(l/iV) contributions. 

When this decomposition is substituted, the self-consistent equation (29) for f^ 2 ' 1 ) is now expressed sym- 
bolically as 



p(2,l) = X 



(45) 



which is diagrammatically represented as 






A 

i i 
i i 
■ i 

>■ S. 



(46) 
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Here, the first and second terms are of 0(1), while the third and fourth terms are of 0(1/N). 

In the present paper, we will not solve Eq. (45) self-consistently, but we will follow the procedure of Rcf. [15] 
in the evaluation of the critical exponent v. That is, we extract contributions logarithmically divergent at 
small momentum, In k, from the LO and NLO terms in r^ 2,1 ) (fc/2, fc/2; fc), which will give the critical exponent 
v when exponentiated. 

Notice that at the LO the equation, T (2 ^ = I + DqT^, is immediately solved by L^ = 1/(1 - D Q ) with 
Do(p) = — (A/6)II(p), which is nothing but the sum of bubble chain diagrams proportional to I(p) in Eq. (34). 
This is easily understood from the 0(1) diagrams shown in Eq. (46). We will denote 1/(1 — Do) with the same 
dashed line as / in the diagrams. From the low momentum behavior of T L <j (fc/2, fc/2; fc) = 1/(1 — Do(k)), 
one should be able to get the critical exponent v at the LO. Indeed, it gives rise to 

T^\kl2,k/2;k) = 1/(1- D (k)) 

~ !*"'<*>= db k (xP- (47 > 

where the polarization LI is evaluated with the scaling form for G(p). There is no Ink dependence in this 
result, but rather it directly gives the exponent v at the LO. Comparing this with the scaling behavior 
fn 2 ' 1 ) ~ £.2-77-1/1/ anc j usm g the LO result for 77, i.e., t^lo = 0, we find that v at the LO is 



= 1 • (48) 



This is the well-known result in the 1/N expansion analysis. Non-trivial correction for v should be obtained 
at the NLO. 

f2 1) 

With r^o = 1/(1 — Dq), we can rewrite Eq. (45) in a form which is more suitable for the perturbation 
expansion: 

p(2,l) = _J_ _J_1 DT (2,1) 

1- D 1 - D N 



1- D V 1 - D J N 
= — h(l+ — - — Do) —Di — - — H , (49) 

where on the third line we have solved the equation iteratively and shown only the LO and NLO contributions 
explicitly. Diagramatically this iterative solution is represented as 






: 1 J + (...) ^ j 

y^K /-\ j — i 

(a) (b) (c) 

The diagrams (b) and (c) correspond to the contributions from 1, and the last two (d) and (e) correspond to 
the contributions from 1 J; D Do in the parenthesis of Eq. (49). We note that \tp 1 operator is first attached to 
(2 1) 

F^o 1/(1 ~ Do) in all the diagrams. 



2. Evaluation of each diagram 

We now evaluate the 0(1/7V) contributions to L^ 2 ' 1 ) as shown in Eq. (50) by using the scaling form (39) 
for G(p). One should note here that the use of the scaling form is a non-pcrturbative prescription. We will 
substitute in G(p) the exponent 77 obtained by the self-consistent KB equation at the NLO . 
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FIG. 2: NLO diagrams (b) and (d) contributing to the exponent v, with explicit momentum assignment. 



Let us examine the four NLO diagrams (b), (c), (d) and (e) in Eq. (50) one by one, seeking for the In A: 
dependence. We note that the common factor T^q 1 attached to \ip 2 operator in all the four diagrams gives 
1/(1 — Do) ~ (6/A^l(?]))fc 1_2r '. The remaining part in each diagram will result in the lnft dependence to 
modify the exponent v. In this subsection we deal with only diagrams (b) and (d) because the contributions 
from other diagrams (c) and (e) are negligibly small as explained in Appendix. 

Firstly, together with the asymptotic form for / (38), diagram (b) is evaluated as 



(b) 



6 



6 



ft 1 " 2 " 
k 1 ^ 



J0^G(p + k/2)G(p-k/2) (^rr^p) 



-2 



d 3 p 



p + k/2\- 2+rl \p - fc/2p 2+ > 1 - 2 M(r ? )- 1 



(51) 



See Fig. 2, for the assignment of each momentum. The minus sign in ^II _1 (p) ~ —I(p) originates from 
that of the second term in Eq. (44). We have introduced an upper cutoff A of the scaling momentum region, 
while we have omitted A v factors, which can be easily restored. Notice that this integral is logarithmically 
divergent when k = 0, which indicates the infra-red dominance in the p-integration and justifies the use of 
the scaling form for G(p). Indeed, the integral yields the lnfc contribution as 



(b) 



6 



XAirj) 
6 

JAW 



fc 1 " 2 " 
fc 1 " 2 " 



d 3 p 



NA( V )J k (2tt) 
1 



3 p 2 +> -2 +> l-2r ; 



ir 2 NA(i]) 



lnfc 



where we have picked up only the most singular part in k — > 0. 

Secondly, diagram (d) (see Fig. 2 for momentum assignment) is evaluated in a similar way as 



(52) 



6 



XA( V ) 
6 

Ufa) 



fc 1 - 2 " 



k 1 - 2 " 2 



NA{iiY 



d 3 pd 3 q 



A 2 {r])Nj (2tt) 6 P 



A/k d 3 ud 3 v 



(2tt) 6 



- 2+r '\p + k\- 2+, i\p + q\- 2+v \p + q + fel-^g 1 " 2 " 



r 2+ ^\u + k\- 2+ ' n \u + w|- 2+ "|« + v + fcr 2+r 'w 1 - 21 ' 



where we have rescaled the variables u = p/\k\, v = q/\k\ and k = k/\k\. The In A; dependence comes from 
two regions of the above integral: (I) \v\ ~ A/ft, \u\ <C |u| and (II) \u\ ~ \v\ ~ A/ft, \u + v\ <C \u\, \v\. 
One notices that the change of the variables, u' = u + v and v' = —v, maps region II to region I and vice 
versa with keeping the integral the same, which means that these two regions give the same In k contributions. 
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Therefore, the In k contribution in the above integral coincides with twice that from region I: 



(d) 



(i 



6 

6 
6 



fc 1 - 2 " 
fc 1 - 2 " 



NA(v) 2 
4 

NA(r]) 2 
4 



2 x 



d 3 u d 3 v 
A/k d 3 ud 3 v 



u- 2+r >\u- 



NA(ti) 

-2 



1T 2 NA(r]) 



2.4(7?) 

lnfc 



(2tt) 6 
K ' k d 3 v 



k\- 2+r >\u + v\- 2+7 >\u 



k\- 2+r >v- 



fc|- 2 +V- 2 " 



(2tt) 3 



(53) 



Finally, we checked that both the coefficients in the diagrams (c) and (e) in Eq. (50) are consistent with 
zero in a numerical integration. Therefore, we do not include these two diagrams as alluded before. Details 
of the calculation are shown in Appendix. Notice that these two diagrams do not contribute to v in the 1PI 
1/N expansion [15]. 



Result 



We collect the NLO corrections (52) and (53) to the LO result (47). As we noticed before, these NLO 
contributions (52) and (53) have the same prefactor r L( j . Therefore, we are able to exponentiate the In A: 
term to obtain ft 2 ' 1 ) at the NLO: 



-(2,1) 
NLO 



fc fc 

2' 2 ; 



:fc 



r 



(2,1) 
LO 



fc fc 

2' 2 

,1-2')- 



fc 



■K 2 NA{rj) 



In k 



XA( V ) 



Then, comparing this with Eq. (15), one finally obtains the exponent v of the 2PI NLO calculation: 



(2PI) 
^NLO 



1 



1 



+ V + ir 2 NA(r)) 



(54) 



(55) 



where i] should be the 2PI critical exponent (43) for consistency. This is our main result. In the limit N — > oo, 
this result of course recovers the LO result (48) together with lini/v_j. 00 r](N) = 0. We plot v as a function of 
N in Fig. 3, where v from the standard 1PI action at the NLO [15] 



v 



(ipi) 

NLO 



1 



3tt 2 
32A 



(56) 



is also shown for comparison. At large N the difference between the IPI and 2PI results diminishes and both 
converge to the LO result t'Lo = 1. At small N the 2PI result stays positive, while the IPI result can become 
negative. We see that the result of the 2PI NLO calculation gives an improved estimate for the exponent v 
at N = 1, 3 than the IPI NLO result. 



V. SUMMARY AND DISCUSSION 



In the present paper we have developed a method to compute the critical exponent v using the 2PI effective 
action. Although v is associated with the diverging behavior of the correlation length £ near the critical point, 
one can compute it on the critical point by analyzing the three-point vertex function r' 2 ' 1 ). Roughly speaking, 
this is possible because r^ 2,1 ) is a derivative of the correlation function with respect to the temperature, and 
thus includes certain information on the deviation from the critical point. In the 2PI formalism we can write 
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1 




FIG. 3: Exponent v from 2PI NLO calculation as a function of N (solid line). The result from 1PI NLO calculation 
(dashed line) and experimental data [22] are shown for comparison. 



down a self-consistent equation for I^ 2 ' 1 ', which is easily derived from the KB equation for the two-point 
function. The explicit form of the equation was obtained to the NLO in the 1/N expansion. We solved this 
equation by iteration to the NLO in the 1/N expansion, and identified from the resultant L*- 2,1 ) the exponent 
v, Eq. (55), as shown in Fig. 3. 

The difference between the 2PI NLO result (55) and the 1PI NLO result (56) comes from two points as 
follows: Firstly, in the 2PI formalism, we deal with the full propagator in contrast to the free propagator in 
the 1PI formalism. Secondly, the sets of the NLO diagrams for r' 2 ' 1 ) are different between the 1PI and 2PI 
formalisms, while there is only one LO diagram which is common in both. Namely, the 1PI NLO calculation 
involves the following five diagrams [15] which should be compared with four 2PI diagrams shown in Eq. (50): 



< 




There is no self-energy insertion in the 2PI diagrams because it is already resummcd in the full propagator. 
In fact, the first 1PI diagram contains one self-energy insertion, and is already included in the 2PI LO 
diagram Eq. (50) (a). This rcsummation of the self-energy diagrams in the 2PI formalism enables us to take 
account of important higher-order contributions into the form of the full propagator, which is the origin of 
the improvement of the 2PI result over the 1PI result in the calculation of the critical exponent. 

Notice that the 1PI NLO result (56) has an apparent flaw at small N. The exponent v must be positive 
v > as it describes the diverging behavior of the correlation length £ near the critical point (cf. Eq. (2)). 
However, the 1PI NLO result becomes negative at small N, although such a small value of N is outside the 
validity region of the 1/N expansion in a strict sense. In contrast, the exponent v in the 2PI NLO result 
remains positive for all N . and it is closer to the experimental data at N = 1,3 [22]. 

Expanding r] of Eq. (55) in N, we see that the 2PI result reproduces Ma's 1PI result (56) [15], and includes 
a part of higher order terms. It shows that the 2PI effective action resums not only leading-log-terms which 
are resummed in the 1PI calculation but also a certain class of the higer-log-terms. 
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q+k/2 

(c) 




FIG. 4: NLO diagrams (c) and (e) 



In the present paper, we did not require the self-consistency for I^ 2 ' 1 ), but rather solved the equation by 
iteration to the NLO. Using conformal invariance in the coordinate space [19, 23, 24] in (28) at the critical 
point, one may analyze the self-consitcnt solution for I^ 2,1 ) to get a better estimate for v. However, one 
should keep in mind that the higher-order calculation of the exponent v in the standard 1PI formalism up to 
C(l/iV 2 ) [18, 20] tends to deviate from the experimental values. Inclusion of higher-order terms by requiring 
the self-consistency for I^ 2 ' 1 ) is an open issue. 
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Appendix 



Here we show details of the calculation of diagrams (c) and (e) in Eq. (50). We closely follow Ref. [15] for 
identification and evaluation of the lnfc contributions. With the momentum assingment shown in Fig. 4, each 
diagram is calculated as follows: 



(c) 
(e) 



6 
6 
6 

JAW 



/fc 1 " 2 " 
fc 1 " 2 " 



(2tt) 3 
d 3 q 
(2tt) 3 



H- l (q)n- L (q + k)G(q + fe/2)T(fe, q) 



n- 1 ( 9 )n- 1 ( 9 + fc)T 2 (fc, 9 ) ( • V x (fc) 



6 A 



d 3 c 



N A(ti) J (2tt) 3 



(58) 



(59) 



where T(fc, q) represents the triangle part and is defined by 



T(k,q) = P//\t + P 

= / ^G(p)G(p + k)G(p-q) 

„ f ^L k - 3+ ^u- 2+ ^\u + k\- 2 +^\u-v\- 2 +^ 
J (27r) J 

= r 3+3 "T (M). (60) 
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Here we have again introduced dimensionless variables u = p/\k\, v = q/\k\, and k = k/\k\. Then two 
diagrams become 



(c) 
(e) 



G 



XA( V ) 
6 

\Adnj 



fc 1 - 2 " 
fc 1 - 2 " 



A? / ^ l - 2 ^ + k\ l -^\v + k/2\-^T {kv) 

4 f d 3 V ! 



NA(rf) J (2tt) 3 



u 1 - 2 "|t> + fe| 1 - 2 "T 2 (fe,u) 



(61) 
(62) 



As we disscuss in the text, we are interested in the lnfc contribution when fc is small. We first note that 
the lnfc contribution appears not from the integration over u in To, but from the integration over v in To. 
Instead, it will appear from the integration over v. However, it is not straightforward to see the power of v 
from the above two expressions. If T)(fc, v) generates V _3 + 3? J 5 then lnfc terms appear in Eq. (61); 



T (k,v) ~ v~ 3+3r ' -> (c) 



A/fc 



dvv 



In fc. 



Similarly, in Eq. (62), if Tq (fc, v) yields v 5+4? ', then lnfc terms appear; 

fA/fc 



Tj(k,v)~v 



-5+4rj 



(e) 



dvv 1 ~ In fc. 



(63) 



(64) 



Below we examine whether these powers indeed appear in T. The definition of To in Eq. (60) can be explicitly 
written as (8 and <p are the angles between two vectors (u, fc) and (u,v), respectively) 

T (k,v) = ^ J u 2 du J dn U - 2+v (u 2 +2ucos6 + l)- 1+n/2 (u 2 + 2uvcos<j) + v 2 )- 1+v/2 . (65) 

To see the power of v in this quantity, it is convenient to perform the Mellin transformation defined by 



(66) 
(67) 



g(s) = / dxf(x)x s \ 
Jo 

fix) = 7T-- / dsx- s g{s). 

Z7T7 J c—iaQ 

The independent part of Eq. (65) can be Mellin-transformed as follows: 



dv(u 2 + 2uv cos (j) + v 2 )- 1+ ^ 2 v 3 - 1 = 2 1 / 2 -"/ 2 (sin 9 i)' ) / 



2-1/ 



1/2 r(s)r(2 - t] - s) /1 + cos 



r(2 - n) 



1 — COS < 



v /4-i/4 



?7 3 7717731 — cos ( 

x 2 Ti -s- - + -,s+ , — - + -, 

2 2' 2 2' 2 2 2 



xu 



= F(s, r), 4>)u s+7J - 2 , (0 < Re[s] < 2 - 77) 

where B is the beta function, is the associated Legendre function and F is the hypergeometric function. 
Here we have used the formula [25] (0 < Re[s] < 2/j, — 7r < 8 < n): 



J {x 2 + 2axcos0 + a 2 )- f *x s - 1 dx = 2^ 1 / 2 (sin(9) 1/2 -'T \p + B(s, 2/i - sJP*/ 2 "^ (cos 0)a 

= 2^- 1 / 2 (si, 



S — 2/i 



-1/2 fl M/ 2 -„ r( s )r(2 M - a ) ^ 1 + cos ^ 1/4 -^ /2 



r(2/x) 



1 — cos ( 



1 1 — cos 8 



X2F1 ( S + ^+-,S-fi+-,fl+-, 



s-2/j, 



(69) 
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Therefore, Mellin transform of Eq. (65) is given as 

/>oo 

T {k,s) = / dvT (k,v)v s 
Jo 



.s-l 



/ u 2 du { dn U - 2+ ^(u 2 + 2ucos9 + l)- 1+ ^ 2 F(s,ri,(f>)u s+ i- 2 



8tt 3 



r^r / dVLF{s,r},4>) I duu s+2 ^- 2 {u 2 + 2u cos 0+1)" 

57T 7 Jo 



87 

If we use the formula (69) again, then Eq. (68) becomes 

' duu s+2r >- 2 {u 2 + 2ucos8 + l)-!+"/2 = 2 1 /2-V2 (sin ^2-1/2 T(s + 2 V ~ !) r ( 3 - 3 '7 - s ) 



(70) 



r(2-rj) 



1 + cos 6* N " ' 



1 — cos ( 



5?7 5 5?7 3 77 3 1 — cos ( 
X2-F1 -s - -7T + r,s+ -r - + 77, 



2 2' 2 2' 2 2' 2 
F'{s,ri,6), (1 - 277 < Re[s] < 3 - 377). (71) 



Therefore, Eq. (65) can be written as 



T o(M = ^ [dSlF(s, V ,<t>)F'(s, V ,e), (l-2r,<Re[ S ]<2-r,). (72) 
Performing the inverse Mellin transformation, we obtain 

T {k,v) = — / dsT (fc,s) V - s 

ZTTi Jc—ioo 
1 i-c+ioo -. 

where 1 — 2?y < c < 2 — 77. For the integration over s, we close the integration path in the right semicircle. 
Then, the poles of F(s, 77, </>) and F'(s, 77, 9) are, respectively, at 

s =2 - T), 3-7?, ■ • • , (74) 
s =3 - 377, 4-37?, • • • • (75) 

Poles we are now interested in are s = 2 — 77 for diagram (c) , and s — 3 — 3?7 for diagram (e) which yield the 
lnfc contributions (see Eqs. (63), (64)). Below we evaluate To(k,v) only at these poles. 
First consider the pole s = 2 — 7/. Since the residue of the gamma function at z = —n is 

f-D" 

lim (z + n)T(z) = ± — f-, (76) 



Residues of F and F' at s = 2 — 77 can be evaluated as 

s - 2 + n)F(s, 77, d>) = 2 1 / 2 - 1 '/ 2 - ( cos ^ ) 



lim (s - 2 + n)F(s, 77, </») = 2 1 / 2 - r '/ 2 - ( cos ^ ' (1 + cos 0) 1 / 2 +"/ 2 = 1, (77) 

->2-?7 2 V 2 J 

lim (,s-2 + 77)F'( S ,77,#) = 2 1 / 2 -"/ 2 (l + cos fl)"/ 2 ' 1 ^ F(1 + ~ 2?/) 

s-v2— »j r(2 — 77) 

Wi^,i+^,iz^y (78) 



17 



Substituting these into Eq. (73), we get 
^3 j dnFis^^F'is^v- 2 ^ 



1 



i r(2 - t?) 

1—3?/ 1 + 3?7 3 — Tj 1— COS#^ -2+r) 



= L lV - 2+ \ 

Next, consider the other pole at s = 3 — 377. Then, 



(79) 



lim (s — 3 + 3r/)F(s, 77, 

5— >2— 77 



lim (s — 3 + 3r))F'(s, 77, 1 



2^(l + cos^-^ r(3 -^ r( ^ + 2??) 



1 . 



5 — 5?; —3 + 577 3 — 77 1 — cos < 



(80) 
(81) 



Therefore, Eq. (73) becomes 



[ ( mF(s,r ]l (j))F'(s,r 1 ,9)v- 3+3ri = [ dcos^ 1 ' 2 ^' 2 ^ + cosfiY 1 

8?r J J 4tH J_ x 



x 2 Fi 



_ 1/2 r(3-3r7)r(-l + 2t7) 
5 — 5?; — 3 + 577 3 — 77 1 — cost 



1X2-77) 

-3+3r; 



EE L 2 7J- 3+31 '. 

As a result, T becomes 

T (k,v) = L lV - 2+v + L 2 v~ 3+3r > + (higer order), 
and the contribution of diagram (c) is 



(82) 
(83) 



(c) 



6 



6 

Jam 

6 

JAM 



fc 1 " 2 " 
fc 1 " 2 " 
fc 1 " 2 " 



d 3 ?j x _ 



2)/ 



N / 8tt 3 



u + fcl 1 - 2 " \v + k/2\- 2+r 'T (k,v) 



4 / >A/fc rf 3 ?; 
^ J ^v^\v + k\^\v + k/2\- 2 ^L 2 v 



-3+3r, 



4L 2 1 A 
A 2tt 2 n ft 



(84) 



Similarly, diagram (e) is 
(e) 



6 



6 
6 

Jam 



ft 1 - 2 " 
ft 1 - 2 " 
ft 1 - 2 " 



NAM J 8tt 3 



piv^lv + k^T^v) 



-4 



J 8tt 3 



-4L!L 2 1 A 
A„4(??) 2tt 2 11 ft 



(85) 



Notice that both diagrams (c) and (e) have the similar structure as those of diagrams (b) and (d). What 
remains is the estimation of the coefficients L\ and L 2 . We evaluated L\ and L 2 numerically, and found that 
L\ ~ 0(1) while L 2 is consistent with zero. Therefore, we conclude that diagrams (c) and (e) could have the 
In k contributions, but are numerically very small, and can be ignored in our calculation. 
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